########## EXERCCIOS 4 - ANA CAROLINA PAVAN #################

### Cervejas ####
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
pref.cervejas <- as.factor(cervejas)
pref.cervejas
class(pref.cervejas)
tab.cervejas<- table(pref.cervejas)
barplot(tab.cervejas,main="preferncia de cerveja por estudantes")
dotchart(as.numeric(tab.cervejas),labels=levels(pref.cervejas),main="preferncia de cerveja por estudantes")
# Barplot tem maior razo dado/tinta!

### Rios ###
data(rivers)
rivers
rios.menores<-rivers<mean(rivers)
rios.menores
tabela.rios<- table(rios.menores)
prop.table(tabela.rios) # a proporo  de ~67%
quantile(rivers,0.75)
mean(rivers)
mean(rivers, trim=0.25)
median(rivers)
#Existe diferena entre os 3 estimadores: a mediana tem o menor valor e a mdia, o maior. A mdia truncada tem um valor bem discrepante da mdia!

### Caixetais ###
caixeta
hist(caixeta$h, main="Distribuio da altura nos caixetais",xlab="altura das arvores",ylab="frequencia")
par(mfrow=c(3,1))
hist(caixeta$h[caixeta$local=="chauas"],main="Altura - Chaus",xlab="altura das arvores",ylab="frequencia")
hist(caixeta$h[caixeta$local=="retiro"],main="Altura - Retiro",xlab="altura das arvores",ylab="frequencia")
hist(caixeta$h[caixeta$local=="jureia"],main="Altura - Juria",xlab="altura das arvores",ylab="frequencia")
## a distribuio de alturas nos caixetais da jureia parece se aproximar de uma distribuio normal; em retiro e chauas, as distribuies de altura so diferentes. Em Chauas tem uma arvore que  "outlier", com 480m de altura, que expande a escala do eixo x at 500. 

### Eucaliptos ###
eucaliptos <- read.table("egrandis.csv", header=TRUE, sep=";")
head(eucaliptos)
class(eucaliptos$rotacao)
eucaliptos$rotacao<- factor(eucaliptos$rotacao, levels=c("1","2"))
class(eucaliptos$rotacao)
library(lattice)
par(mfrow=c(1,1))
bwplot(dap~regiao|rotacao, data=eucaliptos)
qqnorm(eucaliptos$h)
qqline(eucaliptos$h)

### Mais Caixetais ###
head(caixeta)
caixeta$dap = caixeta$cap / base::pi
head(caixeta)
Tabebuia = caixeta[caixeta$especie=="Tabebuia cassinoides",]
plot( h ~ dap, data=Tabebuia ) # assim gero um plot dap/altura para tdos os caixetais
plot(h~dap, data=Tabebuia, subset=local=="chauas")
plot(h~dap, data=Tabebuia, subset=local=="retiro")
plot(h~dap, data=Tabebuia, subset=local=="jureia")
scatter.smooth(caixeta$dap[caixeta$local=="chauas"],caixeta$h[caixeta$local=="chauas"], xlab="DAP", ylab="Altura")
scatter.smooth(caixeta$dap[caixeta$local=="retiro"],caixeta$h[caixeta$local=="retiro"], xlab="DAP", ylab="Altura")
scatter.smooth(caixeta$dap[caixeta$local=="jureia"],caixeta$h[caixeta$local=="jureia"], xlab="DAP", ylab="Altura")
xyplot(dap~h|local, data=Tabebuia)

####### FIM! #####


